Experimental observation of superluminal group velocities in bulk 
two-dimensional photonic bandgap crystals 

D. R. Solli, 1 C. F. McCormick, 1 R. Y. Chiao, 1 and J. M. Hickmann 1 - 2 ]] 

1 Department of Physics, University of California, Berkeley, CA 94720-7300. 

2 Departamento de Fisica, Universidade Federal de Alagoas, 

CO ! CidadeUniversitdria, 57072-970, Maceid, AL, Brazil. 

O 

^ ■ Abstract 

We have experimentally observed superluminal and infinite group velocities in bulk hexagonal two- 
dimensional photonic bandgap crystals with bandgaps in the microwave region. The group velocities depend 
on the polarization of the incident radiation and the air-filling fraction of the crystal. 
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I. INTRODUCTION 



The superluminal propagation of wave packets with faster-than-c, infinite and negative group 
velocities has been observed in a wide range of physical systems, including both passive and ac- 
tive optical media JUS SO, Si- Garrett and McCumber Jq] first predicted that a Gaussian wave 
packet can propagate with negligible distortion through a linear optical medium with a superlumi- 
nal group velocity, if its bandwidth is restricted to a narrow spectral region of anomalous dispersion 
within an absorption or gain line. This prediction was later experimentally verified for the propa- 
gation of picosecond laser pulses through a GaP:N sample O]. Such superluminal behavior does 
not violate relativistic causality. Any analytic function, such as a Gaussian wave packet, contains 
sufficient information in its early tail to reconstruct the entire wave packet with a superluminal 
pulse advancement, but little distortion. Relativistic causality forbids only the front velocity, i.e., 
the velocity of any discontinuity, from exceeding the vacuum speed of light, c. It does not forbid 
the group velocity of a wavepacket from being superluminal. 

The superluminal behavior of electromagnetic pulses can occur not only in passive, absorptive 
optical media, but also in passive, dissipationless optical media, such as periodic dielectric pho- 
tonic bandgap (PBG) structures [ 8] . Optical pulses whose bandwidths are restricted to the PBG are 
exponentially damped not due to absorption, but interference between multiple Bragg reflections 
from planes within the transparent periodic dielectric structure. The superluminal propagation of 
wavepackets through ID periodic dieletric structures has been observed in the optical spectral re- 
gion [9]. Similarly, superluminal propagation has been observed for classical microwave pulses in 
a ID periodic dielectric stack placed inside a microwave waveguide [10]. Other work has demon- 
strated superluminal behavior in the frequency range of tens of MHz in a ID photonic structure 
made of a sequence of coaxial cables with alternating impedance values Jlli . No reported ex- 
periment has yet demonstrated superluminal propagation of electromagnetic wavepackets in 2D 
periodic structures, although we have recently demonstrated for the first time clear experimental 
evidence for the existence of exponential wave decay in such structures at microwave frequencies 
I12I]. Here we report on experimental results which imply that superluminal group velocities do 
occur in such structures. 

Microwave measurements have certain advantages which we exploit here. First, one can per- 
form direct, precise phase measurements using a microwave vector network analyzer (VNA). Sec- 
ond, the wavelength is sufficiently large so that high-precision, well-fabricated 2D photonic crys- 
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FIG. 1 : Experimental setup for microwave transmission and phase measurements. 



tals can be relatively easily constructed. Third, one can more easily tailor the dispersive properties 
of the PBG by changing structural parameters, such as the structure air-filling fraction (AFF). 



II. EXPERIMENT 

The fabrication of our photonic crystals has been described in detail elsewhere Jl^l . In sum- 
mary, we stack acrylic rods in a hexagonal lattice and glue them in place with an acrylic solvent 
cement. We constructed and tested crystals with as many as eighteen layers of rods, and AFFs 
of 0.60, using rods with a 1/2-inch outer diameter. The crystals are "bulk" samples rather than 
"slabs" in the sense that their extruded (non-periodic) dimension is many times the lattice spacing 
in length. 

In Fig. [T]we show a schematic of the experimental apparatus used for the microwave trans- 
mission measurements. We used an HP 8720A VNA, which is capable of measuring both the 
transmission and phase delay of microwaves detected by a polarization- sensitive receiver antenna 
(horn). The VNA was swept from 8 to 14 GHz in 15 MHz increments and was connected to 
identical transmitter and receiver horns. The crystal and receiver horn were placed 1.6 m away 
in a box measuring 60 cm on a side, with microwave-absorbing walls. Microwaves entered the 
box through a 14 cm x 17 cm rectangular aperture. These dimensions were chosen to mini- 
mize diffraction effects through the aperture without allowing leakage around the crystal. With 
the crystal removed and the rectangular aperture closed, the microwave signal at the receiver was 
suppressed by more than 45 dB, indicating good shielding by the box. As previously reported, our 
bulk crystal with AFF equal to 0.60 displays a strong bandgap at « 1 1 GHz ll^ll . 
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FIG. 2: Transmission through an AFF equal to 0.60 crystal with 18 layers, (a) TM polarization; (b) TE 
polarization. 

To establish the reliability of the phase delay measurements, we placed sheets of ordinary 
acrylic of varying widths in front of the receiver horn. The phase delay relative to that of free 
space propagation was linear in the transmitted microwave frequency, with a slope that depended 
on the thickness of the acrylic. Using a simple model of phase delay given by 

A<f>(w) = -d(n(u) - 1), (1) 
c 

where d is the sheet thickness, we were able to extract a value of n = 1.58 for the index of 
refraction of acrylic, in good agreement with the accepted value of 1.61 for this frequency range 

Q. 

The group velocity measurements were carried out by measuring, for both polarizations, the 
spectral dependence of the phase delay for crystals with different numbers of layers. The group 
index was calculated according to the relation [6] 

n g (u ) = n(uo) + — (2) 
where n(u) is the normal (phase) index of refraction and the group index n g is defined as the ratio 
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FIG. 3: Phase index of refraction n(co) for an AFF equal to 0.60 crystal with various layer numbers, (a) TM 
polarization; (b) TE polarization. 

of the speed of light to the group velocity. Thus, superluminal group velocities exist in any spectral 
region where the group index is less than unity. 

In order to calculate the group index, it was necessary to unwrap the phase data obtained from 
the VNA, since it is sensitive to phase modulo 2n. We were able to determine when such phase 
slips occurred by observing that the unwrapped phase delay at fixed frequencies was a monotonic 
function of the number of layers except for sharp changes (phase slips) at certain numbers. We 
adjusted the data by adding 2n x m to each phase delay spectrum, where m is the cumulative 
number of sharp drops observed over smaller layer numbers. We also confirmed the validity of 
this phase adjustment by considering dn/doj, which should go to zero far from the bandgap. We 
found that it did so if and only if we applied this phase delay adjustment as described. 





10 11 12 
Frequency (GHz) 



14 



FIG. 4: Group index of refraction n g {ui) for an AFF equal to 0.60 crystal with various layer numbers, (a) 
TM polarization; (b) TE polarization. 

IH. RESULTS AND DISCUSSION 

In Fig. El we show the measured transmission bandgap of a crystal with eighteen layers and 
AFF equal to 0.60, for both incident microwave polarizations. We consider the plane of incidence 
to be the plane of periodicity (perpendicular to the rods) and define the TE (TM) polarization to 
have its electric field perpendicular (parallel) to this plane. We normalized transmission measure- 
ments by removing the crystal and measuring the total microwave signal at the receiving horn. We 
find that this crystal exhibits anomalous dispersion within the bandgap, passing briefly through 
zero dispersion at the bandgap edges, for both TE and TM polarizations (see Fig. |3j) . We note that 
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the region of anomalous dispersion becomes increasingly well-defined as the number of layers is 
increased. 

From the phase delay data in Fig. |3]we calculate the group index using Eq. El The results 
are displayed in Fig. HJ The region of superluminal n g extends over a wider spectral band for the 
TM polarization than the TE polarization, and deepens more rapidly as the number of layers is 
increased. This result is not surprising given that the TM gap is known to form more rapidly than 
the TE gap as the number of crystal layers is increased II 1211 . 

We did similar measurements (not shown here) using other crystals with AFF equal to 0.32 and 
obtained qualitatively similar bandgaps, although they are narrower and shallower for both polar- 
izations, as expected full , and have slightly different center frequencies. Both crystal sets display 
similar superluminal spectral regions, although we found that those with higher AFF typically 
have greater bandwidth and yield superluminal group velocities which are more superluminal. 



IV. CONCLUSION 



We have experimentally demonstrated that superluminal including infinite group velocities can 
exist for analytic signals whose spectral bandwidth lies within the bandgap of a two-dimensional 
hexagonal photonic crystal. Superluminal phenomena are exhibited by these crystals for both 
tested AFFs, although the effects are more significant (i.e., group velocities become more super- 
luminal and exist over a wider spectral range) for structures of higher AFF. 

This work was supported by ARO grant number DAAD 19-02- 1-0276. We thank the UC 
Berkeley Astronomy Department, in particular Dr. R. Plambeck, for lending us the VNA. JMH 
thanks the support from Instituto do Milenio de Informacao Quantica, CAPES, CNPq, FAPEAL, 
PRONEX-NEON, ANP-CTPETRO. 
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